General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



X-021-75-157 

PREPRINT 


► ' * ' 





POLAR MOTION SPECTRA 
BASED UPON DOPPLER, LP.M,S., 
AND B.I.H. DATA 




MICHAEL A. GRABER 


(N ASA-TH-X“70<»24) POLAB MOTION SPECTPA N75-2797'* 

b/ SED UPON DOPPLER IPMS AND BIH DATA (NASA) 

IS p HC S3.2S CSCL 03B 

UncLds 

(iJ/89 2BB0 7 


JUNE 1975 



GODDARD SPACE FLIGHT CENTER 

GREENBELY, MARYLAND 


X-921-75-157 

Preprint 


POLAR MOTION SPECTRA 

BASED UPON DOPPLER, I. P.M.S. , AND B.I.H. DATA 


Michael A. Graber 
Geodynamics Branch, Code 921 


June 1975 


GODDARD SPACE FLIGHT CENTER 
Greenbolt, Maryland 20771 U. S.A. 


1 



POLAR MOTION SPECTRA 


BASED I’PON DOPPLER, I.P.M.S. , AND B.I.H. DATA 
Michael A, Graber 

ABSTRACT 

Spectra based upon recent Doppler, International Polar 
Motion Service, and Bureau International de I'Heure data are 
compared. An oscillatioii Is found at 1. 3cpy which mif^ht be 
due to an Eulerlan motion of the solid Inner core. An extended 
15-year I.P.M.S. data set is filtered and aiu^lyzed to yield a 
Chandler peak with a period 430.8 solar days and a full width 
at half-maximum of 0.7 days (Q^OO). The data Is reanalyzed 
in overlapping 3-yoar segments and indicates that the exci- 
tation of the Chandler wobble is a dif.crete process and that 
periods as long as three years occur in which the driving 
mechanism is essentially quiescent. 
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I>OL‘\H MOTION SPKCTItA 


HASED i:iON DOPPLER, P.M.S. , AND B.I.H. DATA 

Since Its development In UV17 by Bur^, the maximum entropy method of spec- 
tral analysis (Burg, 19(17; Burg, 1968; Peacock and Treltel, 1969; Lacoss, 1971| 
has pr>' .>d successful in treating short data records which trould not lie resolved 
with the traditional techniques utilizing Fourier transforms or autocorrelation 
functions. The phenomenon of polar motion is one of those which has been suc- 
cessfully treated. International Polar Motion Service (I. P.M.S. ) (Claerbout, 
1969; Currie, 1974] and Bureau International do I'Heure (B.I.H.) (Smylle et al. , 
19731 data sets have been analyzed elsewhere already. Here we develop the 
power spectrum for recent’;/ acquired Doppler satellite tracking data (Anderle, 
19731 and compare it with spectra for B.I.H. and I. P.M.S. data for »hc same 
epoch. In the secoed part of this paper, we develop the power spectrum for a 
filtered segment of I. P.M.S. data. 

The Doppler satellite tracking data spans the four-year period from 1971 
through 1974, with a tw'o-day sampling period, TVo gaps appc*ar in tho data. 

The first is for one data observation for New Year's Eve 1971-1972; the second 
is for 29 observations near the end of 1972. The data set was made continuous 
by constructli\g the appropriate prediction filter and calculating approximate 
values for the missing points consistent with the statistical properties of the 
surrounding data. B.I.H. and !. P.M.S. data were analyzed for comparison 


1 



with thf Doppler data, with data Iuti((th8 of six (1908-1973) and seven (1908-1971) 
years, respectively. 

The maximum entropy technique of spectral analysis was now applied to the 
three sets of data. This technique as described by Burg is utill/.ed by first ge'.- 
eratlng a prediction error filter either by an iterative process which Burg used 
(Burg, 1908; Levinson, 1947| or by a Wiener-Hoph matrix equation. If the 
Wiener-IIoph equation is used, the prediction error filter length is limited to 
10 per cent of the data lengi-h, since the equation uses a matrix composed of 
autocorrelation values which are usually limited to a 10 per cent lag. Then to 
extend the ’Miter further, the iterative technique must be applied, in this Investi- 
gation, It w’as found that application of the iterative process alone gave better 
spectral resolution than the ccmblnatlon of the Wiener-Hoph and Iterative 
-hniques. 

An aspect of the maximum entropy technique, which is not yet well understoml, 
is whethei Linus should be placed on the length of the prediction error filter. 
Given a data set with L points, the filter may be generated to a maximum length 
of L-1. It would seem at first that the maximum length is also the optimum 
length, since the spectral resolution can be estimated as the Nyquist frequency 
divided by some characteristic number of the calculation method, e. g. , the 
number of autocorrelation lags, or Fourier transforms, or elements In the pre- 
diction error filter. Thus, a longer filter would be expected io Improve the 


2 



rosolutior between two oscillations such as the annual and Chardler com- 
ponents of polar motion. 

In practice, as the filter lenRth Is Increased, the spectrum starts to seriously 
deteriorate and often Innirs little resemblance to the anticipated spectral power 
dlsirlbutlon. This point has bc-cn r» co)fnt7a d la^fore [Currie, 1974). Smylle et 
al. (1973) used a 404 point data set, l)ut only (generated a 70 element prediction 
error filter. With 8ri4 data points, Currie 11974) Roncrateo 133 filter elements, 
or a 15 P' 1 ?" cent filter. I'slnK implicit double precision arithmetic and the crite- 
rion descrlbcKll)elow, we have stably extended the filter length to 62-65 per cent 
for three separate data sets having 138, 438, and 725 data points. 

he filter length criterion is based on the fact that each data set contains a 
finite amount of noise power, denoted by H in Fig' •'el. In generating a predic- 
tion error filter, the noise power H is calculated as the noise power estimate P. 
With an infinite data set, as the filler length is increased the « .tlmate Poo would 
slowly converge to the value H, as indicated in the figure ’ descending dashed 
line. However, because of the finite data length L, as the imcr approaches 
length L-1, the noise |K)wer estimate Pj is observed to behave In the pathological 
manner indicated In the figure and <s no longer a good estimate for H, 

In this investigation, the optimum fl’ter length was determined by observing 
the behavior of P^^ and ending the iterative process w'hen the noise power estimate 
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NOISE POWER 



Figure 1. The anomalous behavior of the noise power estimate, Pj^ , 
establlshf s the optimum filter length. For the pole position data 
considered here, filter lengths of 55-65 per cent were utilized. 
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Btartecl to IncreuHv Its rate of c*e8cent. * This criterion has resulte<l in filter 
lengths up to 65 per cent, allowing increased resolution while maintaining spec- 
tral stability. 

The jMmer spectra for the re^d'^n "!• 5 to -t-1.5 cycles per year (cpy) are given 
In Figure 2. One of the more striking features of this figure is the variation In 
resolution of the annual and Chandler )x;aks. In the Oopplcr spectrum, the cen- 
tral region between the two peaks is 32 dB down from the Chandler peak; in the 
I.P.M.S. 8(>ectrum, this value decreases to 25 dB; in the B. I, H. , 17dB. In ad- 
dition the full width at half-maximum for the Chandler peak is 1-1/2 days 

for Doppler, 2-1/2 days for I.P.M.S. , and 7 days tor B. I. H. 

Contrary to expectations, the largest ratio of annual-to-Chandler amplitude 
lo In the Doppler spectrum. It had been thought that the optical data from B. I. H. 
and I.P.M.S. would have been more highly contaminated with annual noise. This 
does not now appear to be the case. The power ratio Is emphasized here -aiher 
than the power estimate for each peak (area under the peak) because Lacoss 
119711 has found that the power ratios are a more accurate estimate in the maxi- 
mum entropy technique. 

The ratio of the period of oscillation to the I’WHM is found to be significantly 
higher for the Chandler oscillation than the typically accepted values of 20-60. 

*in a recent paper Ulrych and Bishop |I97S| discuss a fdter length criterion based upon the mean square 
prediction error. 
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FREQUENCY (CPY) 

Rgure 2. Maximum entropy spectra for threu sets of pole position 
data, showing; the high degree of resolution provided by the Doppler 
satellite data. 
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The vaiut 8 d' termlned from these spectra are: Doppler, 290; I. P. M.S. , 170; 

H. I.H., 00. One consequence of the hl^h ;xjrlod 'KWHM, or Q, for the Chandler 
oscillation is that the energy '> out needed to maintain the Chandler amplitude 
would not be as lar^o as had been previously calculated, Lacoss (19711 has 
shown that the maximujn entropy technique Rives a better estimate of the spec- 
tral shape' than the Blackman-Tukey autocorrelation technique. Thus, a param- 
eter which depends upon the shape, such as Q, shrjuld be better estimated with 
the maximum entropy technlq\ie. 

The spectra In FlRure 2 all have similar structure Ir the rcRlon around 0.0 
cpy. These F>enks are probabij lue to the secular drift of the pole. The B. I. H. 
and I. P. M.S. spectra show a symmetric doublet structure, which indicates a 
linear motion through the pole. A variation in the rate of drift durlnR this jjerlod 
caused the computer program to asslRn a non-zero frequency to the motion. The 
presence of a single strong peak in the Doppler spectrum indicates that the pro- 
gram assigned a finite curvature to the mot.on in ihls short data set. 

All three spectra show a p>eak at approximately -1.0 cpy. This structure, 
noted by Claerbo >t (19691 in I. P. M.S, spectra, indicates a linear component 
or a non-zero elliptlcity In the annual forcing function. It would be caused by 
a linear component in the atmospheric redistribution of mass which accompanies 
the annual seasons. The lack of a similar negative Chandler peak, most apparent 
in the I.P.M.S. spectrum. Indicates that this part of the spectrum is not due to 
the elliptlcity of the earth’s equator. 
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An unexpected peak at approximately 1.3cpy (27<i daya) appears in all three 
spectra. In view at the varyinK data collecting metho^ls associated with the 
Poppler technique ana the optical observations of the D. I.H. and I.P. M.S. , 
and the different sampling periods (2 days. Dopple**; 5 days, U.I.Ii. ; 18 days, 

I. P. M.S. ), this is probably a real oscillation. Tlie period of the peak shows a 
scatter Ir. the three spectra (260, 271, and 283 daya), not considered significant 
since the Chandler and annual periods show a similar, though smaller, scatter. 
This problem of frequency shifts In the application the maximum entropy tech- 
nique to short data records has been noted before (Ulrych et al. , 1973|. 

The power of this 1.3cpy peak is down by two orders of magnitude from 
the Chandler peak in each sp'.ctrum. Indicating that the amplitude of the motion 
is one-tenth of the Chandler, or approximately 1/2 meters. The definition of 
the peak indicated by the period/ I'V'HM is relatively high, being calculated as 
20-40. In the absence of other reasonable explanations for this oscillation, we 
have cons! Jc red the possibility that It Is due to an Eulerian or Chandler-type 
oscillation of the solid inner core. The sense of the circular oscillation is the 
same as the Chandler peak, and thus has the right orientation, in the earth as 
a whole, a 20 km equatorial bulge over a 6400 km radius is sufficient to cause 
a 430 day Chandler motion (from a 300 day Eulerian period). Thus, depending 
on the elasticity of the inner core, all that would be required to generate this 
type of oscillation would be a 4-5 km equatorlsl bulge over the measured 1230 km 
radius of the inner core. Since seismological analyses of the inner core boundary 
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arc uncertntn on the 5 tan level (MiiaMe et al. , 1974), It ia posaib'.c that thla 
bulKo exlnta. The lack of ai^lficant density contrast at the Inner-outer core 
boundary would, however, lead to inefficient coupling to the mantle. 

In view of the excellent .•esults achieved with the I.P, M.8. data for 1968,00- 
1974.85, It wuH decided to Incorporate additional data and reanaly/.e the I. P. M.S. 
polar motion spectra for the 15 year period 1900.00-1974.86. The results of this 
analysis using the maximum entropy technique are shown In Figure 3, labeled 
the unflltered spectrum. As would be expected, the results were not a signifi- 
cant improvement over the previous analysis of aeven years of data. 'Fhe Chandler 
peak is at 437 days with a FVVIIM of 3 days. 

At this point a new technique was utilized which is a variation of the pre- 
dictive filtering methods described by Ulrych et al. 11973) and by Smylle et al. 
11973). The data w’hich Is composed of 298 complex values is predicted in a 
forward and backward direction by a prediction filter of length 160 until 2048 » 2" 
points nave been listed. These points are then analyzed with ihe Sande-l'ukey 
fast Fourier transform algorithm. In the technique of Smylle et al. , the Fourier 
amplitudes of the annual peak were reduced to the background level and an inverse 
fast Fourier transform was performed to obtain a new data set with which to 
generate a filtered power density spectrum. We felt that a better representation 
of the ttower structure of the Chandler peak would be achieved if the Fourier am- 
plitudes ofthe annual peak were multipliefl with the appropriate sine or cosine 
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FREQUENCY (CPY) 

Figure 3. Filtered and unflltered spectra for I. P. M.S. data (1960,00- 
1974.85). The Insert shows the smooth, symmetric structure of the 
Chandler peak with a full width at half maximum of 0. 7 days. 
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fi ’ctlon, and these time series were then subtracted from the original data. 

An)’ Jlscontlnultios in the Chandler motion would be better represented by this 
technique. The resulting cower spectrum ir shown in the lo.ver half of Figure 5. 
Some noise resulting from the subtraction can be observed at 30 dB below the 
Chandler peak, but the addition of double-precision si .w and cosine time series 
with frequencies l.OOScpy and 0. 995cpy can be expected to have very little effect 
at the Chandler frequency 0. 85cp> in the maximum entropy technique. The im- 
pi .^vement observed is due to the removal of the interference effects from two 
nearly eoual amplitude oscillations. The beat period of the annual and Chandler 
oscillaMons is 13 years. Thus wltli 15 years of data the removal of one of 
the oscillations should greatly enhance the resulting spectrum for the remaining 
term. 

As can be seen from Figure 3, the resulting spectrum does show a signifi- 
cant Improvement over the unflltered spectrum. The peak appears at 430. 8 
solar days with a I’WUM of 0. 7 days. This corresponds to a period/ FWHM (=Q) 
=600. The insert in Figure 3 shows the smooth symmetric structure of the 
Chandler perk. 

The value for Q associated with this filtered spectrum is significantly higher 
than commonly accepted values for Q, 20-60. A Q 600 would allow the damping 
of the Chandler motion to be accomplished by the mantle without requiring damp- 
ing by the oceans. According to Stacey [19701, if the Chandler motion is damped 
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by the mantle alone, then the Q of the Chandler oscillation should .^e 7. 5 times 
greater than the Q of the mantle. Q(Chandler) ^600 would Imply (^mantle) -80 
at 0. 85cpy. Analyses of seismic oscillations with frequencies around one cycle 
per minute yield values for the average Q(mantle) between 200-250 (Stacey, 19691 . 

Figure 4 shows the depe.idence of various parameters of the Chandler peak 
for the filtered spectrum of Figure 3 on the length of the prediction error filter 
utilized. The length considered optimum Is indicated with the vertical line. The 
Important points are the general stability of the parameters near the cut-off and 
the change in the behavior of the parameters after the cut-off Is exceeded. 

With the annual oscillation filtered out of the data, adequate rcaolutlon of 
the Chandler oscillation should be possible with the maximum entropy technique 
using approximately two cycles of data. The filtered data was reanalyzed in 
nine overlapping sets of 2.9 years duration. The resulting spectra are shown 
in Figure 5. An interesting result of this computation is the absence of an an- 
nual peaK r. om all nine spectra. The filtering was done using the Fourier am- 
plitudes for the entire 15 year set of data. Variations in the annual oscillation 
would appear as small annual peaks In the nine spectra in Figure 5. The absence 
of any annual peaks Indicates a high degree of stability In the annual forcing 
function. 

The Chandler peak Is seen to vary In frequency and i-’W’IlM from spectrum 
to spectrum. Figure 6 shows the variation In peak period. Note that the average 
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Figure 4. Dependence of Chandler peak parameters on the length of 
the prediction error filter. The parameter values are stable as they 
approach the filter length determined by the error power criterion. 
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Figure 5. Power density spectra for nine overlapping sets of datn. 
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Figure 6. Variation of the peak period and full width at half maximum. 
Note the non-random grouping of the points. 
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of ti'ese nine Chandler periods Is 433.6 solar days, which Is very close tc 430.8 
days calculated from the filtered spectrum In PM^ure 3. This Indicates that these 
peaks are all associated with the Chandler oscillation but are perturbed by some 
statistical fluctuation centered on the Chandler frequency. Figure 7 shows the 
variation of FWHM with the shift of the p«iak fronr. 430.8 day a. Seven of the points 
In this figure seem to show a clear functional dependence of the FWHM on the 
distance from 430.8. 

To determine the extent of this functional dependence, the least squares 
criterion was used to generate a quadratic line of regression. As was expected, 
when all nine points are used the correlation Is roor. When the two wild points 
are dropped ((shlft=33. 6, FWHM=2. 5), (19.0 25.1)), the following line of re- 
gression Is established 

Y = FWHM (days) 

X = period shift (days), greater than zero 

Y = 0. 711 + 0. 129 X + 0. 010 X^ 

Correlation Coeff. = 0. 97 

Note that this line of regression estimates the FWHM at zero period shift (430.8 
days) to be 0.71 days. The FWHM calculated from Figure 3 Is 0.7 day.>. The re- 
gression analysis recovers the FV/HM exactly to the accuracy of the Initial value. 

Recall that the mean of the nine peaks, 433.6 days, recovered the Chandler 
period to less than 1 per cent. If the two values dropped In the regression analysis 
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Figure 7. A possible function.il dcpondencfe between the full width at hell’ 
maximum and the period shift. Seven of ihe nine points establish a quad- 
ratic line of regression with a correlation coefficient of 0. 97. The two 
omitted points are based on two adjacent overlapping data seta. 


17 



are removeH from the mean, the new mean is 432.3 days which misses the 15 
year Chandler period b” 3.5 parts in a thousand. 

We ijclieve that there is a strong possibility that the variations in peak fre- 
quency and I’ M shown in Fijjure 5 are due ti small discontinuous phase shifts 
in the Chandler »^sclllatlon. These phase shifts are caused by changes in the 
moment-of-inertla tensor of the earth (earthquakes, etc.) which result in a 
change in the location of the figure axis. They would change the Chandler oscil- 
lation from Figure Sa to either Figure Hb or 8c. it can be seen that a phase shift 
greater than zero in Sb would cause a spectral analysis program transforming a 
short (lata set to assign the spectral power to a higher frequency than is in fact 
occurring In the oscillation; and in 8c, a phase shift less than zero would cause 
the analysis program to compute a peak frequency lower than the actual frequency. 
In both cases with non-zero phase shifts, the break in the time series would cause 
a broadening of the peak. 'Hils association of frequency shift and line broadening 
is the type oi behavior shown in Figure 7. 

If this pustulated phase-shift process takes place very frequently (f > 1 cpy) 
and in a random manner, then the correlation observed in F’igure 7 for the 2.9 
year analysis would not exist. I'he multiple random phase shifts would cause 
a period shift of zero over the 2.9 year analysis. In the 15 year analysis this 
randomness does nullify the period shift. 'Fhe high correlation observed in 
Figure 7 indicates that in the seven short data sets used in the regression analysis 
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Figure 8. A possible mechanism for the period shifts observed in the spectra, 
(a) X-component of unperturbed Chandler motion, (b) and (c) A shift in the 
figure axis could result In the X-component undergoing this type of motion. 
Comn;ter analysis would yield period shifts and line broadening. 





probably only one siKnificant phaan shift is occurring in each period. In the two 
wild periods, 1964. 35-19(17. 20 and 1965-80-1968.65, probably more than one 
phase shift has occurred. And it is possible that during the periods 1968. 70- 
1971.55 and 1971.60-1974. 45 the driving mechanism was quiescent or the phase 
shifts occurred near the ends ct the period. This second possibility is very likely 
for at least one of these two periods since the overlapping period 1970. 15-1973.00 
exhibits a significant shift and broadening. 

This general type of behavior in which the driving force for the wobble is 
applied discretely and randomly once every two or three years was postulated 
by Dahlen (1971 1. 

An attempt was made to correlate the FWHM with numbers of large earth- 
quakes using lists provided by the National Earthquake Information Service. No 
correlation was found. Probably the absolute number of earthquakes is too simple 
a quantity to correlate with the I'WHM of the Chandler oscillation. 
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FfGURE CAPTIONS 


Fltfure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 
Figure 6. 


Figure 7. 


The nnomaloue t)etuivior uf the noise power estimate, , establishes 
Uie optimum filter length. For the pole position data considered here, 
filter lengths of 55-65 per cent were utilized. 

Maximum entropy spectra for three sets of pole position data, showing 
the high degree of resolution provided by the Doppler satellite data. 

Filtered and unflltorcd spectra for I.P. M.S. data (1960.00-1974.85), 
The Insert shows the smooth, symmetric structure of the Chandler 
peak with a full width at half maximum of 0. 7 days. 

Dependence of Chandler peak parameters on the length of the predic- 
tion error filter. ITie parameter values are stable as they approach 
the filter length determined by the error power criterion. 

Pow'er density spectra for nine overlapping sets of data. 

Variation of the peak period and lull width at half maximum. Note 
the non-random grouping of the points. 

A possible functional dependence between the full width at half maxi- 
mum and the period shift. Seven of the nine points establish a quad- 
ratic line (H regression with a correlation coefficient of 0.97. "nie 
two omitted points are based on two adjacent overlapping data seta. 



KiKure 8. 


A possible mechanistii for the period shifts observed In the spectra, 
(at X-component of unperturbed Chandler motion, (b) and (o) A shift 
in the figure axis could result In the X-component undergoing this 
type of motion. Computer analysis would yield period shifts and 


line broadening. 


